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Abstract. Criteria for the membership of individual meteors in meteor streams are 
discussed from the point of view of their mathematical and also physical properties. 



^ ' Discussion is also devoted to the motivation. It is shown that standardly used criteria have 

on : 

CD ■ unusual mathematical properties in the sense of a term "distance" , physical motivation 
■ and realization for the purpose of obtaining their final form is not natural and correct, 

o : 

CD and, moreover, they lead also to at least surprising astrophysical results. A new criterion 

O : 

"--■^ i for the membership in meteor streams is suggested. It is based on probability theory. 
Oh' Finally, a problem of meteor orbit determination for known parent body is discussed. 
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5t 1- Introduction 

Meteor streams are composed of meteors originating from a parent body (comet, asteroid) 
during an ejection process. It is very important to know the properties of a meteor stream 
also in order to make some possible conclusions about physical properties of the parent 
body, ejection process(es) and time of possible ejection(s). One of the most fundamental 
feature of a meteor stream is the membership of individual meteors, i. e., if an individual 
meteor belongs to the given meteor stream or if it is a background meteor. Meteor stream 
membership criteria were suggested for this purpose. 

Fundamental characteristics of a body in the Solar System is its orbit characterized 
by orbital elements. This holds also for meteors. Orbital elements of meteors originating 
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from one parent body are similar. Thus, it is natural that the most simple meteor stream 
membership criteria are based on investigation of meteor orbits and their orbital similar- 
ity. Some other criteria may be also used, if it is possible: e. g., physical composition of 
individual meteors, based on the observed spectra. 

D-criterion of Southworth and Hawkins (1963) is the most frequent meteor mem- 
bership criterion used in literature. It measures the orbital similarity of two individual 
meteors and thus it enables to find complete meteor streams. Any meteor is represented 
by a point in a phase space of orbital elements. The phase space is a "five-dimensional or- 
thogonal coordinate system" (in this section quotation marks are used for the text taken 
from Southworth and Hawkins, 1963) and each element is considered as a coordinate. 
Southworth and Hawkins take the following orbital elements: q (perihelion distance), e 
(eccentricity) , i (inclination) , ft (longitude of the ascending node) , u (argument of peri- 
helion) . According to the authors, "the distance between two points is a natural measure 
of the difference between the two corresponding orbits" . The orbital similarity function 
- "distance between two points" (meteors) A and B - D(A, B) is then defined as 

[D(A,B)} 2 = (q A - q B f + (e A - e B f + (<l sin 
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(1) 



where I AB is the angle made by the orbital planes given by the formula 



„ . Iab\ 2 ( ' . ia - «s\ 2 , . . . . { n A - n B \ 2 

2 sin— — 1 =12 sin I + smi A smi B 12 sin I , (2) 

and tt ab is the difference between the longitudes of perihelion measured from the 
point of the orbits intersection 

o • / f i A + i B \ . f n A - n B \ I AB \ 

ir AB = uj a - lo b + 2 arcsin <^ cos I I sin I I sec > , (3) 

or, to a sufficient approximation, when i A and i B are small, 
n AB = (Q A + u A ) - {Q B + u B ) . (4) 

This D-criterion is used for identification of meteor streams: if D(A, B) < D c , where 
D c is a constant assumed as a threshold value, then orbits of meteors A and B are similar 
and both meteors may be members of the same meteor stream. 

A little modified D-criteria were suggested by Drummond (1981) and Jopck (1993). 

The aim of this article is to analyze mathematical and physical properties of 
Southworth-Hawkins' D-criterion (properties of Drummond's and Jopek's D-criteria are 
analogous). Finally, we present a new meteor stream membership criterion. 
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2. Mathematics and D-criterion 

In the construction of their D-criterion, Southworth and Hawkins were originally inspired 
by the well-known definition of the distance in Euclidean space. As a consequence they 
have obtained D-criterion in the form (1) (and also Eqs. (2)-(4)). As it is considered, D- 
criterion measures the distance between orbits of two meteoroids (meteors). However, if 
it is so, then D-criterion (1) must fulfill properties required for a quantity called distance. 
The standard properties of a distance are closely connected with the so-called metric 
space. Definition states: 

Let X be a set with elements u,v,w,.... A nonnegative function p defined on the 
Cartesian product X x X is called a metric if it satisfies the following axioms: 

(i) p(u, v) = if and only if u = v ; 

(ii) p(u,v) = p(v,u) ; 

(iii) p(u,v) < p(u,w) + p(w,v) . 

A set X with a metric p is called a metric space. 

(Metric - distance. The property (iii) is called the triangle inequality.) 

Now, question is, if these properties are fulfilled also for D-criterion (1). One can 
easily verify that triangle inequality is violated. It means that triangle inequality, which 
is an evident property of a distance, is not fulfilled in the case of measuring "distances" 
between meteor (meteoroid) orbits. 

As an evident property of the D-criterion we introduce the following one. If D(u,v) 
is smaller than D(u,w), then orbits u and v are more similar than the orbits u and w. 
However, due to the violation of the triangle inequality, the orbits v and to may be more 
similar than one would expect on the basis of general conception about distance: 

< D(v,w) < D(u,w) - D(u,v) . 

We can formulate this in terms of meteor orbits: the distance between meteors u and v 
is small, the distance between meteors u and w is large, but the distance between meteors 
v and to may be small. 

From the mathematical point of view it would be useful to have a D-criterion which 
fulfills the properties (i)-(iii). It may have also other forms than presented by Eq. (1), i. 
e., it may contain not only terms of the form (aj — ctk) 2 , but also, e. g., terms of the 
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form \ctj — a>k\, and so on. We know from mathematics that the number of possible (from 
a mathematical point of view) metrics is infinity: if p(x,y) is a metric, then pi(x,y) = 
p(x, y)/ {1 + p(x, y)} is also a metric. One of the most simple modifications of Eq. (1) 
satisfying properties (i)-(iii) may be obtained by substituting the term (e^ + es) /2 by 
any constant; the consequence of the modification (e^ + es) /2 — > 1 is, e. g., that the 
group of Cylids (Southworth and Hawkins, 1963, p. 271) consists of two meteor streams 
{Dsn = 0.06, 0.10, 0.10, 0.10, 0.14; D new = 0.24, 0.19, 0.59, 0.59, 0.77). 

As for meteors, standard procedure is to choose one orbit as a reference orbit. The 
reference orbit may be given by mean values of the known members of a given meteor 
stream. Any meteor A is a member of the stream if D(A, M) < D c , where M represents 
an orbit defined by mean values of orbital elements and D c is a constant assumed to be 
a threshold value. Since D(A, M) is given by Eqs. (1) - (4), we may have D(A, M) < D c 
and D(B, M) < D c for some two meteors A and B of the stream, but D(A, B) > 2 D c . In 
other words: the first two inequalities assert that meteors A and B have similar orbits but 
the last inequality states that their orbits are not similar! (Points A and B are situated 
inside a sphere with a centre M and a radius r, but the distance between A and B may 
be greater than 2r!) This is the consequence of the triangle inequality violation. (If D — > 
and, moreover, Ae/e -C 1, sin Att/2 <C 1 (i. e., orbits are identical or almost identical), 
then the violation of the triangle inequality plays no role. However, this may not be the 
case occuring in astronomical applications (see Klacka and Volosin 1996).) 

Those, who are interested in mathematical properties of semi-metric (or, even in a 
more general case of semi-pseudometric) defined by Eq. (1), we refer to section 18 in Cech 
(1966). 

3. Physics and D-criterion 

Southworth and Hawkins (1963) present also physical arguments for the choice of D- 
criterion in the form of Eqs. (l)-(4). However, their arguments are not convincing and 
thus one should take Eq. (1) as an empirical criterion. We will justify this statement now. 

The physical model of Southworth and Hawkins is based on the idea that the change of 
metcoroid's orbital elements with respect to those of the parent body may be represented 
as an average value of the changes (weighted by velocity) during one period of the parent 
body. However, this idea has no advantage. Nor it is a simple idea (hypothesis), nor does 
correspond to real processes: the meteoroid is ejected at once (at one moment) from the 
parent body. 
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Moreover, mathematical realization of the physical model is incorrect and some other 
physical inconsistencies arise in the process of mathematical calculations. According to 
the idea, we should write for the change of any orbital clement G 

AG = (V„G) • (At;) . 

However, the authors use in their calculations 
AG = |V„G| \Av\ , 

which is different from the correct value, since, in general, Av is not in the direction of 
V V G: 

-1 < {(V„G) • (A«)}/{|V„G| \Av\} < +1 . 

The consequence of these three equations is that even if the model of the authors 
would be correct, its realization is incorrect - they add together not changes of orbital 
elements, but their extremal values with positive signs at each instant. Although this 
would seem acceptable in the sense of calculations of maximum possible changes - which 
is not the author's interpretation, and, even if it would be -, it is incorrect: for the same 
meteoroid, Avr may be dominant for one orbital clement, Av$ may be dominant for 
another orbital element (see Eqs. (A8) in Southworth and Hawkins (1963); we used the 
same notation as Southworth and Hawkins). The other important nonphysical result rests 
in the assumption that \Av\ is proportional to circular velocity at the distance r: 

\Av\ oc U , i.e., \Av\ oc r -1 / 2 . 

(In the appendix of Southworth and Hawkins, there should be 
AG cx J r- 1 ' 2 |V„G| dt . 

Fortunately, calculations of the authors correspond to this result.) 

I am not sure about the correctness of the final results of Southworth and Hawkins 
(their appendix and Fig. 1), e. g., my result is: lim e _o {(l/ a ) Grad q} = 0.64, which 
seems to be not consistent with results depicted on Fig. 1 in Southworth and Hawkins' 
article. (The mean values calculated in appendix of the authors are calculated from 
absolute values of quantities, so as the mathematical lemma could be used.) 

According to Fig. 1 in Southworth and Hawkins (1963), the authors make a conclusion 
that their D-criterion may be applied for e < 0.85. In reality, it is used even for e ~ 1. 
Individual terms of the sum in Eq. (1) are comparable in their values also for e > 0.85, 
which is not consistent with Fig. 1 in Southworth and Hawkins. This also shows that Eq. 
(1) is not consistent with the physics suggested by the authors. 
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4. Astronomy and D-criteria 

On the basis of the previous two sections we can conclude that meteor stream member- 
ship criteria (Southworth and Hawkins 1963, Drummond 1981, Jopek 1993) have two 
unpleasant properties, as for mathematics and physics: 

i) triangle inequality does not hold for arbitrary orbits, 

ii) physics of the criteria is unknown. 

If we want to make a simple physical model, we may imagine that mctcoroids are 
ejected only at parent body's perihelion and calculate the change of meteoroid's orbital 
elements with respect to those of the parent body. One can obtain 

Av \ 2 _ (p A - Pm) 2 {&a - e M ) 2 



vpbJ 8 P 2 m 8 (1 + e M ) 

(smI AM f + [ - CM — ) (simT AM f , (5) 
\1 + eM J 

where vpb is parent body's speed at perihelion (ejection of the meteoroid A), subscript 
M corresponds to the parent body; p = a(l — e), a - semimajor axis. We may define, 
on the basis of Eq. (5) the following quantity: 

[D(A,M)f = a (p A - p M f + P (e A - e M f + 

(&mI A M) 2 + 7 (sin tt am) 2 , (6) 

where a, (3 and 7 are numerical constants for a given meteor stream: a = 1/(8 p M ), 
P = 1/(8 (1 + ejvf) 2 ), 7 = e M /(l + cm) 2 - The first two terms of the sum are equal in 
the problem of two bodies. In reality however, they are often not comparable (some terms 
of the sum are often negligible in comparison with the others) for real meteor streams. 
This shows that although physics of Eq. (6) is simple, it does not correspond to useful 
meteor stream membership criterion. 

Practical advantage of the Southworth and Hawkins' criterion in comparison with the 
criterion defined in Eq. (6) is that the individual terms of the sum are more comparable 
than it is in the case of Eq. (6). However, if this is the only requirement for the choice 
of the criterion for practical usage, then we can write better criterion at once (numerical 
factor £ = cm): 

D(A,B) = \q A - q B \ + \e A - e B \ + 2 ] sin %A ~ %B \ + 2 £ | sin * A ~ WB \ . (7) 

In practice, the quantity eM should be calculated as the mean value of eccentricities 
of bodies forming given meteor stream. Individual terms of the sum of Eq. (7) are more 
comparable than it is in Eq. (1) (moreover, triangle inequality is also fulfilled; remark: 
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Eq. (6) does not represent any distance - the index M cannot be chaanged into B !). 
Many modifications of Eq. (7) may be used, e. g., \qa — Qb\ — » \qa — <7b|/(<Za + the 
same substitution for eccentricities, their various combinations, etc.. 

There is no reason for the fact that individual terms of the sum in Eqs. (1) or (7) 
should be statistically comparable for a given meteor stream. On the contrary, one should 
expect that some terms are dominant for some meteor streams, other terms may be more 
important for other meteor streams. Terms of the sum should be weighted. 

There are other unpleasant properties of the membership criteria of Southworth and 
Hawkins (1963), Drummond (1981), Jopek (1993). All of them use perihelion distance q 
as an orbital element. The consequence is that meteors with large dispersions in semi- 
major axis a may be classified as members of a given stream, which does not seem to be 
real - there is very small probability that one or several meteoroids can be ejected with 
velocities much greater (even in orders of magnitude) than the other meteoroids from 
the same parent body (or bodies of comparable properties). But, when we would like to 
make a substitution q — > a, or, p — > a (e. g., in Eqs. (1) and (7)) the individual terms in 
D(A, B) are much less comparable and usually the term with a is dominant (probably, 
weighting coefficients might solve this situation). The last remark: many people using 
D-criteria incorrectly calculate the mean values of angular quantities. 

In conclusion of this section we can make a statement that it is not possible to give a 
simple D-criterion based on physical arguments. Thus, one can write many forms of D- 
criteria. If we would apply various D-criteria on meteor streams (and their background) 
containing several tens of meteors, we should expect more than 80 % (or, perhaps, even 
more than 90 %) coincidence between them. If we would compare them with the results 
obtained by Tisserand parameter, the coincidence may even decrease to about 50 - 70 
%. (Of course, it has no sense to use various D-criteria for meteor streams containing less 
than 10 known members - various D-criteria lead to completely different results.) 

All statements made in the last paragraph should be understand in the way that it 
is not very useful to use D-criteria based only on the problem of two bodies. Tisserand 
parameter, invariant of the motion in the restricted three-body problem, may be a better 
criterion. However, one must bear in mind that meteoroids are perturbed on their orbits 
not only by gravitational forces. Nongravitational forces may also be very important. And 
some of them may be of stochastic nature. This is the reason of our suggestion presented 
in the following section. 
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5. Probability, Statistical Mathematics and D-criteria 

Since wc do not know any simple physics which can define in a simple way a meteor 
stream, we take the data (set of orbital elements for various meteors) as a random sample. 

Since four orbital elements completely define a meteor stream (intersecting the orbit 
of the Earth), we will take four orbital elements Qi, Q2, Q3, Qi- Let the real distribution 
of orbital elements of meteors in the stream may be approximated by density function / 
(Qi, Q2, Q3, Qi) = f(X). If we define the meteor stream as a set of bodies with led, 



P(X e ft) = / f(X)dX = a , 



(8) 



then there is a probability less than 1 - a that objects with X e Q' belong to the 
stream. 

6. Meteor Orbit Determination 

Another problem concerning the orbits of meteors is the problem of finding radiants for 
given parent bodies. In determining orbital elements for a meteor corresponding to the 
meteoroid A initially ejected from some parent body B, we must use the following sets 
of equations: 

dF d ( p A \ n . , , _ PA , , n . 



d[3iA d(3iA i 1 + e A cos u A J 1 + e A cos lu a 

corresponding to the case if the ascending node intersects the orbit of the Earth, and 

9F ' aJ-|- ^ }=0,i = lto5, 2± = 1) (10) 



dfiiA dfliA 1 1 — eA cos oja J ' 1 — cos lua 

for the case if the descending node intersects the Earth's orbit. j3iA correspond to 
orbital elements of the meteoroid (body A), the constraint corresponds to the fact that 
meteoroid may strike the Earth (circular orbit of the Earth is supposed, for the sake of 
simplicity) . The quantity F may be of the type 

F = a x (E A - E B f + a 2 (L A - L B f + 

a 3 (v T A ~ v TB ) 2 + a 4 (v RA - v RB ) 2 + 

a 5 (v T fA - v TfB ) 2 , (11) 

where 

2 _ / 1 c a 1 C B 



(Ea — E B ) = 

PA PB 

(La - L B ) 2 = pa + Pb - 2 \fp~AP~B x 



{cos^a — n B ) siniA sini B + cosi^ cosi B } , 
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(v TA - v TB f 



(1 + e A f (1 + e B f 

Pa Pb 
2 (l + e A )(l + e B ) xX 



y/PA PB 

X = cos(Qb — &a) {sinujA suilub + 
COSUJA cosujb cosi A cosz_b} + 
sin(Os — Q A ){smoj A cos ujb cos is — 
sin lob cosuoa cosia} + 
cosuoa cos u b siniA sin is , 
(vua v RB f = e l+ e l - 2 ^LxF, 

PA PB ^PA PB 

Y = cos(fis — n A ) {cosuja cosuob + 
sh\to A sinus cosi A cos is} + 
sin(£ls — £Ia){svclu) A cos ujb cosz^ — 
smu)B coscj^ cosis} + 
smoj A sin lob sini^ sin is , 

(vrfA - vtjb) 2 = A - — V 

PA 

(1 - e B /V2) 2 _ 

PB 

(1 - e A I V2) (1 - e B I yg) 
, x ^ ? 

z = cos(ri B — o^) x 

{(sinwA — cosw^) (sinews — cost^s) + 
(sinwA + cqsuja) (sinews + cosujb) x 
cosi A cosz_b} + 
sin(£!s — Q A ) x 

{(sinwA — cosuja) (siu^b + cosujb) x 
cos is — 

(sinwA + cosw^) (sinews — cost^s) x 
cos 2.4} + 

(sinwA + cosw^) (sinus + cosujb) x 

sini^ sinzs . (12) 

One of the coefficients ot\ — 0:5 may be put equal to 1 (e. g., a\ = 1). a— coefficients 
are functions of orbital elements and their values must be determined from known pairs 
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"parent body - meteor stream" . As for the physical sense of the individual terms of the 
sum in Eq. (11), they correspond to: energy, angular momentum, perihelion velocity and 
radial velocity of maximum value, transversal velocity for true anomaly / = 3 ir / 4 in the 
problem of two bodies. The a-coefficients guarantee that Eqs. (9)-(12) can be applied on 
real situations. The advantage of the form of Eq. (11) is that it contains known physical 
quantities. In principle, we can choose other forms. However, the requirement that they 
must contain all five orbital elements in an independent way may not suffice in obtaining 
good coincidence between theoretical and observed radiants. Local and global minima 
may be important, in general. 

7. Conclusion 

We have shown that standardly used method for determining meteor stream membership 
is not correct from the point of view of mathematics, physics and astronomy. We have 
presented correct method, based on probability theory and statistical mathematics. 
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